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Abstract 

We study inclusive semi-leptonic B — > X c £i/g decay using the power counting 
iTic ~ -y/AQCD^b- Assuming this scaling for the charm-quark mass, the decay 
kinematics can be chosen to access the shape-function region even in b — ► c tran- 
sitions. To apply effective field theory methods in this region we extend SCET to 
describe massive collinear quarks. We calculate the tree-level decay rate, including 
0(AQci)/mb) power corrections, and show that it factorizes into a convolution of 
jet and shape functions. We identify a certain kinematical variable whose decay 
spectrum is proportional to the universal leading-order shape function familiar 
from b — > u decay, and speculate as to whether information about this shape 
function can be extracted from data on b — » c decay. 



1 Introduction 



The standard tool in heavy quark physics, the expansion in inverse powers of the heavy 
quark mass, has been extended recently to describe the kinematical situation in which 
the hadronic jet produced by the decaying heavy quark carries a large energy of order 
mb, but has a small invariant mass of order ^m b AQ CI ). Although processes occurring 
in this so-called shape-function region of phase space do not permit a local operator 
product expansion, they can be described by an expansion in terms of non-local operators 
evaluated on the light-cone [1, 2]. This expansion is most conveniently discussed within 
the framework of the soft-collinear effective theory (SCET) [3, 4, 5, 6], as has been done 
in studies of soft-collinear factorization in b — » u transitions up to one-loop order [7, 8, 9] 
and including l/m& power corrections [10, 11, 12, 13]. 

Our goal in this work is to examine the effects of final-state quark masses on inclusive 
semi-leptonic B decay occurring in the shape-function region. For the up, down, and 
strange quarks, these mass effects are presumably very small. This is clearly not the 
case for the charm quark, which is in fact most often treated as a heavy quark with mass 
m c ~ nib. in this work we suggest a different power counting for the charm-quark mass, 
namely m c ~ a/Aqcd^- This power counting is satisfied numerically and allows us to 
access the shape- function region even in b — > c transitions [14, 15]. For decay kinematics 
restricted to the shape-function region the charm quark can be treated as a collinear 
quark with finite mass. It is a simple task to include such a collinear quark mass in 
the SCET Lagrangian [16, 17], but this has not yet been done beyond leading order or 
applied to a phenomenological example. 

In this paper we formulate a version of SCET appropriate for describing processes 
involving massive collinear quark fields interacting with soft and collinear gluons, and 
use it as a tool to calculate the hadronic tensor for inclusive semi-leptonic b c decay 
in the shape-function region. We calculate the hadronic tensor at tree-level, including 
0(AQCT>/ m b) power corrections, and show that it factorizes into a convolution of jet and 
shape functions. The structure of this factorization formula mirrors closely that of the 
b — > u case, the differences being that a slightly different set of shape functions is needed 
at sub-leading order and that the shape functions depend on a different combination 
of kinematical variables as compared with the massless case. We identify a certain 
kinematical variable whose differential decay rate is directly proportional to the universal 
leading order shape function, and speculate as to whether information on this function 
can be extracted from a study of b — > c decay. This variable is a generalization of the 
P + variable suggested in [18], which has been reconsidered in the framework of SCET in 
[9,19]. 

The organization of the paper is as follows. In the next section we discuss our 
power counting and outline the calculation. We construct the SCET Lagrangian for a 
massive collinear quark in Section 3 and derive the corresponding transition currents in 
Section 4. In Section 5 we use this version of SCET to calculate the hadronic tensor 
for semi-leptonic b — > c transitions in terms of a factorization formula. Differential 
decay distributions are discussed in Section 6 and compared with the well-known local 
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expansion of inclusive b — > c spectra in Section 7. In Section 8 we discuss possible 
phenomenological applications of our results and we conclude in Section 9. 

2 Power Counting and Factorization in B — ► X c £p£ 

We are interested in calculating decay distributions for inclusive B — > X c £z/^ decay in the 
shape-function region using effective field theory methods. In this section we will define 
our power counting and explain how a factorization formula arises at tree level. The 
central quantity of interest is the hadronic tensor, which contains all the QCD effects in 
the semi-leptonic decay and can be used to derive the differential decay spectrum. We 
define the hadronic tensor as 

W» v = -lm(B(v)\T^\B(v)), (1) 

7T 

where we use the state normalization (B(v)\B(v)) = 1. The correlator T^ u is given by 

T^ = i J d 4 xe- iq - x T{J^(x)J u (0)}. (2) 

Here J M is a flavor- changing weak transition current, which is 

j» = cY(l- l5 )b (3) 

for B — > X c £u£ decay. We shall use the notation of [20] and write the hadronic tensor in 
terms of five scalar structure functions 

= Wi {VpVv + v„p v - g^vp - ie^a p a v p ) (4) 
- W 2 + W 3 v^v v + W 4 (p^v u + v^p v ) + W b p^p u , 

where the independent vectors are chosen to be v, the velocity of the B meson, and 
p = m b v — q with m b the b quark pole mass and q the momentum of the outgoing lepton 
pair. We use the convention e 0123 = — 1. 

The hadronic tensor is expressed in the effective theory as a double series in the 
perturbative coupling constant and a small parameter A, which we define through the 
relations 

Aonn . o m 2 



^E^A 2 and ^f~A 2 . (5) 
m b ml 

This choice correlates the two small scales m c /m b ~ a/Aqcd/^6- Quantum fluctuations 
are described by three widely separated scales, m\ ^> m 2 A 2 ^> m 2 A 4 , provided that the 
jet momentum p satisfies p 2 — m 2 ~ m 2 A 2 . We will refer to these scales as the hard, 
jet, and soft scale respectively. We expect that the fluctuations occurring at these scales 
can be factorized into a convolution of hard, jet and soft (shape) functions, similarly to 
the massless case [7, 21]. We will outline the steps used to prove such a factorization 



2 



formula to all orders in perturbation theory, under the assumption that there are no 
additional complications compared to the massless case, but do not give a complete 
proof for decay into charm quarks. An essential part of a factorization proof would be to 
show that the IR divergences associated with the renormalization of the leading-power 
jet function are the same for massive and massless SCET, and that the IR divergences 
of the new operators at sub-leading power can be treated in dimensional regularization. 
Our aim here is to explain how our explicit tree-level results can be interpreted as the 
leading-order term of such a presumed all-orders factorization formula. 

We calculate the hadronic tensor by means of a two-step matching procedure, QCD — > 
SCET — > HQET. In the first step, we remove fluctuations at the hard scale by matching 
the QCD Lagrangian and transition current onto their corresponding expressions in 
SCET. We derive the SCET Lagrangian in Section 3 and discuss the transition currents 
in Section 4, to relative order A 2 in the power counting. In addition to these SCET 
expressions, we also need the sub-leading HQET Lagrangian 



C 



(2) 

HQET 



1 



2mb 



h v {iD s fh v + g^Ka^F^K 



(6) 



In general, C mag ^ 1 due to fluctuations at the hard scale, and the SCET transition 
currents are multiplied by Wilson coefficients depending on quantities at the hard scale. 
These Wilson coefficients define the hard functions in the factorization formula. However, 
the SCET transition currents and HQET Lagrangian have trivial hard coefficients at tree 
level, so we will not need to discuss these hard functions any further. 

After this first step of matching, the correlator (2) is calculated using the effective- 
theory expressions for the Lagrangian and currents. This will be the subject of Section 
5. The effective theory contains a collinear charm-quark field and an HQET field for the 
b quark, as well as soft and collinear light quarks and gluons. In the next section we 
will define more precisely what we mean by collinear and soft, for now we simply note 
that the collinear fields fluctuate at the scale m 2 A 2 , whereas the soft and HQET fields 
fluctuate at the scale m 2 A 4 . Since the B meson contains only soft degrees of freedom, 
the matrix element of each term in the correlator factorizes into the form 

{B\Tf\B) = (B | h v [soft fields]^ | B) <g> (0 1 [collinear fields] 1 0) 

= Si <8> Ji , (7) 

where the <g> stands for a convolution. The vacuum matrix element of the collinear fields 
defines a set of perturbatively calculable jet functions Jj containing fluctuations at the jet 
scale. The tree-level jet functions are given in terms of propagators of the collinear fields. 
Calculating these jet functions removes the collinear degrees of freedom and defines the 
second step of matching SCET — > HQET. The matrix element of the soft fields between 
the B meson states is calculated in HQET and defines a set of non-perturbative shape 
functions S^. The differential decay distributions are then written in terms of the Wj, 
each of which factorizes into a convolution of jet and shape functions. 
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3 Matching onto the SCET Lagrangian 



In this section we derive the SCET Lagrangian for soft and collinear gluons interacting 
with collinear charm-quark fields carrying a mass m c ~ Am&. The leading-order result 
has been given previously in [17]. Here we extend the derivation to sub-leading order, 
using the position space formalism of [6]. We do not discuss the Yang-Mills Lagrangian, 
because its matching onto the effective theory is the same as in [22]. 
The starting point is the QCD Lagrangian for massive quark fields, 



where iD^ = id 11 + gA^. The derivation of the effective theory proceeds in the usual 
way. We introduce effective theory fields according to the momentum scaling of their 
Fourier components. In terms of two light-like vectors n + ,n_ satisfying n + n_ = 2, the 
components of a collinear momentum p^ scale as (n + p c ,p c ±,n_p c ) ~ m b (l,A, A 2 ) and 
those of a soft momentum as p^ as ~ mb(A 2 , A 2 , A 2 ). 

Following the standard procedure, we decompose the QCD gluon field into a collinear 
field A 1 ^ and a soft contribution A%. With our power counting for the jet momentum, 
we only have to consider charm-quark momenta which satisfy p 2 — m 2 ~ m 2 A 2 . These 
are described by massive collinear quark fields. Because of the finite charm-quark mass, 
loop diagrams do not contain soft divergences related to charm propagators. Therefore 
charm quarks appear neither in the soft Lagrangian nor in the 5-meson states. This 
is in contrast to the b — > u decay, where the effective Lagrangian contains interactions 
between soft and collinear light-quark fields. 

Notice that in the region of phase space where p 2 c — m 2 ~ m c Aq CD ~ m 2 A 3 , we 
would have to consider another effective theory which describes exclusive hadronic charm 
resonances. Since A/m c ~ A is a small quantity, one can switch to HQET and split 
the charm-quark momentum into p c = m c v' + k', where k' is a residual momentum. 
Because the energy transfer to the charm quark is still large, it is boosted by a factor 
of 0(1/ X) with respect to the 5-meson rest frame. Therefore, the charm velocity scales 
like v' ~ (1/A, 1,A), and the residual momentum scales as k' ~ ra&(A, A 2 , A 3 ). One can 
then introduce an effective theory field h v > carrying a dynamical momentum k', which is 
described in terms of a "boosted HQET" . Restricting ourselves to collinear quark fields, 
we have to make sure that the considered phenomenological observables are not sensitive 
to the exclusive region, see the discussion in Section 8. 

For the collinear quark fields we define the projections 



£qcd = ip(iI/>-m c )ip, 



(8) 



4 



^c(x), 



7)(x) 



4 



i) c (x). 



(9) 



Inserting these into (8) and solving the field equation for 77, one finds 



V = 



4+ 1 



(iip L - m c )£, 



(10) 



2 in+D 
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so that the Lagrangian becomes 

i 



£qcd = i 



in-D + (ilp± - m c )- -(ilp± + m c ) 

in + D 



2 t (id 



The manipulations performed so far are just a rewriting of the QCD Lagrangian. In 
the next step we expand this Lagrangian as a series in the parameter A. The expansion 
is related to the soft gluon fields contained in the covariant derivative iD^ = id^ 1 + 
gA^ + gA%, as well as the multi-pole expansion of these soft fields. The result for the 
leading-order Lagrangian is 

Cf = £ (in_D + (il? ±c - m c )^—(iip ±c + m c )) ^±£, (12) 

while the sub-set of power-suppressed terms proportional to the mass is 



4- = m2 ^^r +As ^D^ 

4- = m &i^ e ^i^/^Tt> 

4Ss = mMx±d±9As), (13) 

where £ < - 1 - ) and contribute to relative order A and A 2 to the effective action. In 
the terms above the collinear fields are evaluated at x, but the soft fields are multi-pole 
expanded and depend only on x^_ = (n + x/2)nt = x + n^L. We will not make this explicit 
in the following, where it is always understood that 

h v (x) = h v (x + n"_) , A s (x) = A 9 (x+nt) . (14) 

The remaining power-suppressed terms are the same as in the massless case. One can 
still make the standard field redefinitions (£ = Y^°\ etc.) to decouple the soft gluon field 
from the leading-order collinear Lagrangian (12), a fact which is crucial to factorization. 

The Lagrangian can be written in a manifestly gauge-invariant form by redefining 
the collinear fields using the methods developed in [22]. Applying these techniques, we 
find that the result for the leading-order Lagrangian (12) is unchanged, but to relative 
order A 2 the power suppressed terms are the same as in the massless case and read [22] 



4 1} =l x ^_W c gF; v W)^i 

4? = \t(n-x)nW_W e gF^W!^t, 

4? = \^± P n-W c [D p ±s ,gF^W^, 
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+ i^^^f^^'^^^xc^e (15) 

The power-suppressed mass terms first appear at relative order A 3 . To see this, we note 
that after making the field redefinitions and working in the collinear light-cone gauge 
n + A c = 0, one effectively replaces 

1 1 



in a- D in+d 



0(X 6 ) (16) 



in the expansion of (11) (see eq. (23) of [22]). It is then easy to show that all mass terms 
not contained in the leading-order Lagrangian are suppressed by at least a factor A 3 . 

While no explicit power-suppressed mass effects appear to C(A 2 ) accuracy in the 
effective Lagrangian after making the field redefinitions, we will see later on that they re- 
appear in the time-ordered products of SCET Lagrangians and currents. The calculation 
with either form of the Lagrangian and currents yields the same result, namely the 
appearance of one new sub-leading shape function as compared to the massless case. 



4 Effective theory currents 

In this section we consider the matching of the flavor-changing weak transition current 
defined in QCD onto its corresponding expression in SCET. For b — > c transitions, this 
matching takes the form 

(WOqcd - z- mbVX [J^ + + J m + •••] , (it) 

where ip c is a collinear charm-quark field. In general, the Jj are convolutions of a short 
distance Wilson coefficient with operators built out of SCET and HQET fields. In what 
follows, however, we will work only at tree level, where the convolutions reduce to simple 
multiplication. The power-suppressed currents we need here can be deduced from the 
results in [22] by noting that at tree level all powers of the mass are related to the n 
field. Using (10) in conjunction with [22], we find for the mass terms 

•ft = rn e & l —W c Th v , (18) 

1 —in + D c 

42 = l —W c Tx ± D ±s h v , (19) 

1 —in + D c 

42 = m c ^^^r^[iip ±c Wc}h v . (20) 
z in + D c Arn h 

The convention is such that derivatives do not act outside of the square brackets, and 
the soft gluon field and the heavy quark field are multi-pole expanded. We have again 



6 



used the field redefinitions of [22] to write the result in a manifestly gauge invariant form. 
The remaining terms are the same as in the massless case [22, 23, 24]. The sub-set of 
these that will be needed later on when calculating the tree-level time-ordered products 
is 

J(°) = lW c Th v , 

= tw c rx ± ^ ±s h v , j 2 (1) = -ti^—^w&rhv, 

in + D c 1 

J? = lW c T^- n + D s h v , 4 2) = lW c Y^^D^ s D» ±s h v , 

jf = -liq ±c —^wJ±Tx^D^ s h v , Jf = £W C T ^h v . (21) 
in+D c 1 Zm b 

5 The hadronic tensor at tree level 

In this section we obtain an expression for the hadronic tensor in terms of a factorization 
formula. We work at tree level and include 0(\ 2 ) power corrections. We start by 
calculating the time-ordered products appearing in the correlator (2) using the effective 
theory Lagrangian and currents derived in the previous sections, setting A c = as 
appropriate at tree level. 

Many time-ordered products are possible at 0{\ 2 ). It is convenient to divide these 
into two groups: those that vanish when m c — > 0, and those that do not. Since the 
sub- leading Lagrangian (15) in its gauge-invariant form is actually the same as in the 
massless case, we can identify the time-ordered products which do not vanish when 
m c — > from eq. (54) of [12]. For these terms the difference between the two cases lies 
only in the leading-order Lagrangian, and we can account for this difference by replacing 
the propagator for massless collinear quark fields by that for massive fields, which is 



(0\T{Z(x) aa Z(y) bf} }\0) = tA(x-y)5 ab 




a/3 



f d 4 k c _ ik ( x _ y ) in + k 
J (2vr) 4 k 2 -m 2 c + ie 

The Latin (Greek) indices refer to color (Dirac) indices. The function A(z) satisfies 

in-dA(z) = 5 w (z) + — *— (ml - (td ± ) 2 ) A(z). (23) 

in + o 

We will always work in the frame where p± = 0, such that in tree-level expressions d± 
can be dropped. 

At tree level the factorization formula contains up to a double convolution between 
the jet and shape functions. The tree- level jet functions are related to (products of) 
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X 





(a) 



(b) 



Figure 1: Representative Feynman diagrams involving soft gluon emission. These give 
rise to (a) bi-local and (b) tri-local shape functions. 



propagators, and some of the convolution integrals can be performed explicitly, after 
which the results can be written as a convolution over a single variable only. To show how 
this works and to shorten intermediate expressions, we introduce the integral operators 



[m c = m c — it, n + p > 0) 



I 2 * f = — / d xe 



(2tt)< 



dn_k 

< 

2tt 



k 2 — ml' 



-in-kx-L 



POD 

= i / dx + e i{n - p - m " /n+p)x+ f(x+), 
Jo 



n^k — m1/n + p 



/(*+) 



(24) 



- j dx + dz + e m - px + J 



dn-k dn-k' ^_ in _ k ( x+ _ z+) ^_ in _ k , z+ 
2tt 2vr 

1 \ f I 



n + p \ri-k — m? c /n + p J \ri-k' — m? c /n + p 



h * — — 

n + p 



j dz + g(x + ,z + ). (25) 
Jo 



In terms of the factorization formula, the functions / and g represent generic shape 
functions depending on one or two variables respectively, and the momentum-dependent 
propagators correspond to tree-level jet functions. These are represented by Feynman 
diagrams such as Figure 1(a) for f(x+), and Figure 1(b) for g(x+,z+). Performing the 
integral to arrive at the final lines of (24, 25) is equivalent to performing a convolu- 
tion integral in the factorization formula. In (24) this is trivial, but in (25) the object 
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f* + dz + g(x+, z+) introduces an "effective shape function" of a single variable only. To 
keep the structure of factorization as presented in (7) clear, we will be careful to dis- 
tinguish these effective shape functions from the shape functions defined directly by the 
multi-local matrix elements of soft operators between B meson states. 

The integral operators I2, h are the same as in [12] after replacing ri-p — > ri-p — 
m 2 c /n + p in the exponential of J| defined there. We can thus obtain the part of the 
correlator that does not vanish when m c — > from eq. (54) of that work by making 
this replacement, and leaving out the term involving the soft quark field, which does 
not appear in the massive case. We also need to identify a small number of terms 
which vanish when m c — > and cannot be derived from previous results. First, we need 
insertions of the transition currents J m . Using the short-hand notation 

j\j B iC c = i j d 4 xe ipx T ^f A (x)J B (0)i J d A zC c {z)^ , (26) 

and working in soft light-cone gauge n_A s = 0, the relevant time-ordered products are 

■J (0)t 4? + J^J {0) = —h * Wrx(o), (27) 

n + p 

WW = * K(x)r ^r^(o), (28) 

J^J (0) = 0. (29) 

(Recall that the soft fields are multi-pole expanded, see (14).) The only additional mass 
terms are related to insertions of 



£?> = -iin-x) n^_W c gF s ^W^, jj 2) = ^^-n + D s h v . (30) 
They contribute to the correlator through the time-ordered products 



j(2)t j(0) + j(o)tj(0) i£ (2) = -_^£-/ 2 * (ix + )h v (x)in + "dr 7 ^rh v (0) 

* h v (x) gn + A s {z)T^Th v {0) 

n + p 2 

im 2 c f x+ - <— -rjb. 



i} yh* j dz + h v (x)(-in + D s )(z)T—Th v (0). 



(31) 



We used an integration by parts to write the field strength tensor in terms of the soft 
gluon field and then applied (23) to simplify (31). The simple form of the tree-level 
propagators has allowed us to replace n^x with n + xml/(n + p) 2 in this equation, as is 
easily verified by writing the factors of x in momentum space as derivatives acting on 
tree-level propagators. These expressions can be made gauge invariant by inserting the 
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appropriate Wilson lines. After including these Wilson lines, the tree-level result for the 
correlator can be written entirely in terms of multi-local covariant derivatives of the type 
in (31). In an arbitrary gauge these are defined by 



(YUD»Y)(z) ab = id»5 ab + (Yi[iD»Y])(z) ab . 



(32) 



The indices refer to color and the derivative does not act outside the square brackets. 
At tree level the result for the correlator with arbitrary Dirac structure is 



T^iu) = - J 



dx+dn_k 



,i(u— n_fc)x_)_ 



(h v Y)(x)F 
+T 



2tc 
4- 



n-k + ie 



x 



+ 



+ 



m 



2 4. 



T(Y^h v )(0) 



2 ' n + p (n + p 2 ) 2 
{Ch v Y)(x)T r ^T(Y^h v )(0) i j d*z£ { mET (z)} 

+ ^[(h v (-i^ s )Y)(x)r 7 ^r(Y^h v )(o) + {KY){x)r^r{Yiiip a h v )(oj\ 

(^(-^^)r)(x)f^ 7M± ^r(rt^)( ) 

+ (W(x)f^ 7 ^r(yt ( _^Mx ) ^ )(0) ^ 



n + p 



+ 



n + p 



- r ( iz + (^y)(x)rt(rt ( _ l ^ i)H ^ l) y )( , )r( yt^ )( o ) 

fP Jo 2 



+ 



im c 

( n+P y 



~dz + {h v Y)(x)(Y\-in + D s )Y)(z)T^T{Y^h v m\, 



(33) 



where we have defined u = n_p — m 2 c /n + p. 

In order to calculate the hadronic tensor, we need to take the imaginary part of the 
correlator and decompose the matrix elements between B meson states in terms of scalar 
functions. The matrix elements of the operators on the final two lines of (33) could be 
used to directly define a set of effective shape functions depending on the variable u. To 
keep factorization explicit we do not take this approach, and introduce effective shape 
functions only in terms of convolutions involving multi-local shape functions defined 
through B meson matrix elements. After contracting with the color indices of the tree- 
level jet functions, we need the following matrix elements (the Greek subscripts indicate 



Dirac indices and e 



ie 



lit/ pa 



n p _v <T ) 



(B\Ch v Y)(xUY^h v )(0)p\B) = ^ n 



S(x+), 



(34) 



pa 



(B\(h v Y)(x) a (Y j h v )(0) l3 i I d*z C% ET (z)\B) 



1 1 (I +i 



2m b 2 



s(x + ), 



(3a 
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(B\(h v Y)(x) a (YHD^Y)(z)(Y^h v )(0) p \B) = 

i (i + i 



(35) 



2 V 2 



0a 



iS'(x+)v» + [iS'(x+)-T 1 (x + ,0)+T 1 (x + ,z + ))n l L 



pa 



t(x + )-f 2 (x + ,0)+f 2 (x + ,z + )\ , 



(s|(^y)(x) a (yt<£>' ± <i>' x y)(«)(yt/i w )(o) /J |B> 
& ( i + y 



(36) 



/3a 



^ 3 (a;+,^+)- 



/3a 



The results are written in momentum space using the convention 

S,(x 1+ ,...,x n+ )= [ du 1 ...du n e- i ^ Xl +"-^ x ^S i (u 1 ,...,u n ) (37) 



for the Fourier transform of a generic shape function Si. We now introduce effective 
shape functions defined by the convolutions 

u s (u) = / du 1 duj 2 J 2 (u;u 1 ,u 2 )[u 1 (cu 1 )5(uj 2 )+U 1 (uj 1 ,uj 2 )], 



ti(u) — dcuiJ 2 (u;ui,0) 



cuiS(coi) — 2 J du'Ti(ui,u') 



+2 J duj 1 duj 2 J 2 (u;u 1 ,cu 2 )T 1 (uj 1 ,uj 2 ); 



(38) 



J 2 (u; oji, u 2 ) = Im 



{n+pf 



7T (p^ -m 2 c )(pl l2 -m 2 c ) 
= —(S(u uj 2 ) - 8{u - ivi)), 

U! 2 



(39) 



where Pu>i2...j = P — ^i2...j^+/2 and oui 2 ...j = uj\ + oj 2 + . . .ujj. J 2 is the tree-level jet 
function related to the product of propagators in (25). 

The contributions to the hadronic tensor which do not vanish when m c — > can be 
obtained from eq. (59) of [12] by replacing n_p — > u. The terms which vanish when 
m r — > are 



Will = I— ^.r 
n + p 2 

ml 1 



i + i 



IT 



+ 



mi 1 



-tr 



(n +P ) 2 2 



-tr 



2 2 



M 2 2 



i±if^±r 

2 2 



(40) 
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To proceed further, we use that the Dirac structures for the semi-leptonic decay in the 
Standard Model are 

f = 7^(1-75), r = 7"(l- 75 ). (41) 

The correction linear in m c vanishes for this V-A structure. Performing the traces, and 
including now all terms, we find that the components of the hadronic tensor are given 
by 

2 



W 1 



n+p 



[(- 




- *) 




n+p J 



S(u) 



S(U 



2m b 



uS(u) — t(u) u s (u) + u a {u) m\ 



m b 



n + p 



(n+p) 



Wo 



Wa = 



2u 



(n+p)' 



n+p 
t(u), 



S(u), 



(n+pf 



uS(u) — t(u) t(u) + u a (u 



m b 



+ 



n+p 



(42) 



The appearance of the effective shape function t\(u) within the tree approximation is 
specific to decay into charm quarks, as is seen by setting m c — > 0. This shows that the 
effects of the final-state charm-quark mass are not merely kinematical, but introduce 
new non-perturbative structure into the factorization formula. This new structure is 
relevant only in the shape-function region. In Section 7 we will extrapolate our results 
to the OPE region by performing a moment expansion, and find that the moments of 
t\(u) are determined entirely by those of the leading order shape function S(u), so that 
in this region the effects of the quark mass are purely kinematical. 



6 Differential Decay Spectra 

From the scalar components Wi of the hadronic tensor one can express the triply differ- 
ential decay rate as [20] 

1 d 3 T 



12r d(n_p)d(n+p)dx 



(43) 



(n+p — n-p) { (1 + x — n^p — n+p)(n+p + n~p — x — n+pn_p)—^-W\ 



+ (1 — n+p + n+pn-p) — W2 + [x(n+p — x) — n+pn-p] — (W3 + 2mbW 4 + m b W^) >, 

where r = G 2 F ml\V c b\ 2 / I92ir 3 and x — 1 — 2Ei/m b . Here and below hatted quantities 
are normalized to m b , i.e. p = p/m b , m c = m c /m b . We have used that only W\ contains 
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a leading order term in A to expand the differential decay rate to 0(\ 2 ). The integration 
range is given by 

* 2 
771 

— c - < n^p < x < n + p < 1. (44) 
n + p 

Any decay distribution can be obtained from (43) along with the result for the hadronic 
tensor in (42). Of special interest for semi-leptonic b — > u transitions is the singly 
differential spectrum in the variable P + = n_p + (M B — m b ), because of the possibilities 
to extract |V^&|[19]. An examination of (42) shows that u = n_p — m 2 /n + p is the natural 
partonic analog of P + for b — > c decay , since this is what appears in the argument of 
the shape functions. (We leave the discussion of an appropriate hadronic variable to 
Section 8.) Changing variables and performing the integrals over n + p and x, we find 

1 dT / 14 m _ m 2 \ „, . s(u) . m 2 



— — — = 1 8^ S(u) + -^-A— UAu) 

T du \ 3 m b mi J 2m b m^ 

+-^—[t(u)+u a (u)-5u s (u)}, (45) 
Smt 

where the allowed phase space is < u < m b — m 2 c /m b . This result agrees with [11] in 
the limit m c — > 0. 

A comment is in order concerning the limits of integration used in calculating (45). 
After changing variables, the integration region is given by 

u H < x < n + p, — + - v M + ™c < n +P < 1- (46) 

Tlj^p Z A 

The collinear expansion is valid in the region n + p ~ 0(1), and breaks down for the lower 
limit above, which corresponds to n + p ~ m c ~ 0(A). However, the doubly differential 
rate behaves like A 2 in this region of phase space. Hence the contribution to the u 
spectrum from this region is in total of order A 3 and does not affect our analysis. 

We can also derive the lepton energy spectrum in the endpoint region, where an 
experimental cut limits the lepton-energy variable x to values satisfying x ~ 0{\ 2 ). 
Such a cut restricts the other two kinematic variables to the shape-function region, 
since m 2 c /n + p < n-p < x and m 2 c j[x — u) < n + p < 1 ensure that n^p ~ 0(X 2 ) and 
n + p ~ 0(1). Integrating over n + p, we find that the leading-order result for the doubly 
differential spectrum is 



1 d 2 Y 



2S{u) 



3mt 2m b r 



[x — u) 2 (x — u) 3 



+ 0{\ 2 ). (47) 



Tq dxdu 

We have also calculated the 0{\ 2 ) power-suppressed terms but do not quote them here. 



7 Comparison with the local expansion 

A more familiar treatment of inclusive b — > c decay uses the operator product expansion 
(OPE) to calculate the total rate [25]. Within this approach, the singly differential spec- 
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trum in the variable u at leading order in a s and including power-suppressed corrections 
up to dimension 5 is given by 

-— = (l-8— 

T du \ ml 



We have not written terms proportional to \i i2 5(u), since these are A 2 jm\ corrections 
which cannot be recovered from our calculation in the shape-function region, which was 
carried out only to 0(A/mi l , m^/ml). The OPE result is valid in the kinematical region 
where all momentum components are of order nib- We can compare the results calculated 
in the OPE region with those derived in the shape-function region by considering u ~ mj, 
in which case moments of the effective shape functions can be expressed in terms of local 
HQET matrix elements [1] (see also [10, 11, 26]). This expansion is given by 

S(u) = 5(u) - y <$"(«), s(u) = - (Ax + 3A 2 ) 8'(u), 
2A 

u s(u) = — t^S'(u), u a (u) = \ 2 5'(u), 

t(u) = -\ 2 6'(u), tl (u) = -^5"(u). (49) 

The higher moments correspond to operators of higher than dimension 5 in the OPE 
and are neglected. 

We will derive the moment expansion for t\(u) as an explicit example, since it is 
unique to decay into charm quarks. Examining the definition (38), we see that it is 
sufficient to derive the moments of the function Ti(ui, 002), since those of S (00) are known. 
At dimension 5, we only need to consider up to the first moment with respect to oo\ or oo 2 . 
Moments of the form J duiidui2 7i(<^i, ^2) ^1 give a vanishing contribution to ti{u), as is 
seen by inserting Ti(cji,cj 2 ) ~ 5^(001)5(002) into (38). This means that we only have to 
calculate the first moment with respect to oo 2 . To do this, we start with its definition in 
terms of a tri- local matrix element (see (32) and (35)) 

' ^(B\(h v Y)(x)[Y^D;Y](z)(Y^h v )(0)\B) = / doo 1 doo 2 e-^ x ^ 2Z+ T 1 (cu 1 , oo 2 ), (50) 



2 

where the derivative does not act outside the square bracket. The first moment with 
respect to oo 2 is obtained by acting on both sides of (50) with —in^d z and then setting 
x and z to zero to obtain a local HQET matrix element, which is calculated using 

1+f 



{BlKT^iDWKlB} = -tr , . 



2 



(g^ -vW)— + ia> w — 
3 2 



(51) 



We then find J duJiduj2T 1 (uj 1 , oo 2 )oo 2 = 0, as can be seen by using the identity 

-in-d z \Y*iD,y\(z) = n»_[YUgF^Y}(z), (52) 
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and noting that the local HQET matrix element vanishes. We conclude that Ti(ui, u> 2 ) is 
irrelevant for the moment expansion of ti(u). Inserting the moment expansion for S(u) 
into (38), we find t\(u) = —^-5"(u) + . . . , as quoted above. 

Replacing the shape functions in (45) by the moments in (49), we find that the 
local expansion of the shape- function region result reproduces the OPE result (48). We 
performed the same comparison for the lepton energy spectrum (47), including power- 
suppressed terms, and confirmed that the moment expansion reproduces the OPE result 
[25], expanded appropriately for m 2 c ~ x ~ A 2 . 



8 Leading-Order Shape-Function Effects 

Neglecting sub- leading shape functions and perturbative corrections, the dT/du spectrum 
is directly proportional to the leading-order shape function S(u). This raises the prospect 
of extracting information on the shape function from data on inclusive charm decays, a 
question which we will now address. Before performing any phenomenological studies 
we will slightly modify our treatment of m 2 corrections. In (45), we kept only terms up 
to m 2 , as is consistent with the SCET power counting. However, higher-order kinematic 
corrections involving the charm-quark mass arise from phase-space integrals and are 
easily included. Keeping these additional kinematic corrections, but neglecting all other 
sub-leading effects, the spectrum becomes 

~ = f(P)^h (53) 

where /(p) is the usual phase-space function 

/(p) = (1 - 8p + 8p 3 - p 4 - 12p 2 lnp) , (54) 

and we have defined p = m 2 /m 2 . The numerical effect of the phase-space function is 
sizable. In particular, using /(p) ~ (1 — 8p) would differ at p ~ 0.1 significantly from 
the value obtained from (54). However, once the logarithmic contribution is included, 
the numerical result is almost exact. 

In an experimental analysis one has to rewrite (53) in terms of hadronic variables. 
To this end, we propose to study the spectrum in the variable 

U = n.P-^ (55) 

where n±P = n±p + {Mb — writ,) = n±p + A are the hadronic light-cone variables. Since U 
depends on n + P one must make this replacement before integrating over n + P, in which 
case the integration region is 



U 1 /— ~ Ml 

2 + 2~V + D ~ n+P ~ Ms ' ° " U ~ Mb ~ 



Mt 



(56) 
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One has to be careful when writing the spectrum in terms of U because 

rr * M D K 

u = U - A + 



n + P n + P — A 

&U-A+ ® — = + ... (57) 

n + P n + Pn + P V ; 

The last two terms are of order A 3 and A 4 respectively as long as n + P ~ M B , but become 
order A 2 at the lower limit of the n + P integration, when n + P ~ M B . However, one can 
show that the contribution of this region to the U spectrum is power-suppressed, so we 
can ignore this subtlety and write u = U — A. The shape function is then a function of 
S(U — A), and the decay distribution is simply 

l§ = fip) S(U-A). (58) 

At leading order p = ml /ml xs Mp/M B . Eq. (58) opens the possibility to at least cross- 
check the results for the shape functions obtained from heavy-to-light decays. However, 
a few comments are in order. 

Unlike in the case of heavy-to-light semi-leptonic decays, almost 80% of the total 
inclusive rate is already exhausted by the two exclusive decays B — > Dti>i and B — > 
D*£u£. This may indicate that the real world is not very close to the limit A — > using 
the power counting we are suggesting. For the exclusive decay B — > D£u e the U spectrum 
is concentrated at U — 

W^ -nB-DWHU) (59) 

while for the exclusive decay B — > D*iv^ the decay rate is a function of 

= Ml, - Ml ^ X*Ml 
n + P y n + P 

Even when n + P ~ Mp ~ AM^ we have U ~ X 3 M B and this is beyond the sensitivity of 
the leading-order approach, since the smearing due to the shape function is of the order 
A 2 . Thus the contribution of both ground states may be written as 

mB^m^ + = [r{S ^ D(p() + r(S ^ D . ep()] m (60) 

In terms of a moment expansion of the decay rate, we see that these two exclusive decay 
modes contribute only to the zeroth moment, and conclude that all higher moments of 
the U spectrum can originate only from excited or non-resonant states. For this reason, 
it would be interesting to see if and how well the relations for the moments, such as 

S ivv k% =Knp) - /^§=B A ' +A ') /W (61) 

are satisfied. In particular, we find that to leading order in the SCET expansion the ratio 
of moments should be the same as the corresponding ratio of P + moments in B — > X u £ue 
or of photon energy moments in B — > X s / -f. 
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9 Conclusions 



In this work we examined inclusive semi-leptonic b c decay using the power counting 
m c ~ a/Aqcd^ for the charm-quark mass. With this power counting the decay kine- 
matics can be chosen to access the shape-function region even in the decay B — > X c £p£. 
We applied effective field theory methods by modifying SCET to include a massive 
collinear charm quark. This defines a consistent power-counting scheme in terms of the 
small parameter A ~ m c /mb ~ a/Aqcd/^6- We matched the Lagrangian as well as the 
weak transition current at tree level and including sub-leading terms up to 0(\ 2 ), and 
used these to derive the hadronic tensor in the shape-function region. 

The results are similar to those for b — > u transitions. This led us to identify a certain 
variable which is analogous to the variable P+ in b — > u decay. As in the heavy-to-light 
case, at leading power and at tree level the singly differential spectrum in this variable is 
directly proportional to the leading-order shape function. We showed how the moment 
expansion of this differential spectrum matches the local OPE result in an appropriate 
limit. 

It remains to be checked whether the power counting we have applied for the charm 
mass yields a consistent picture or whether the conventional counting m c ~ rrib is more 
appropriate. A crosscheck can be performed by measuring the moments of the distribu- 
tion in this variable and by comparing them to those obtained from the P + spectrum 
in B — > X u £u£ and from the E y spectrum in B — > X s ^. Our leading-order prediction is 
that these moments should be equal up to the phase-space factor which appears in the 
total partonic rates. 
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